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ABSTRACT 

It has been argued that a substantial fraction of massive stars may end their lives 
without an optically bright supernova (SN), but rather collapse to form a black hole. 
Such an event would not be detected by current SN surveys, which are focused on 
finding bright transients. Kochanek et al. (2008) proposed a novel survey for such 
events, using repeated observations of nearby galaxies to search for the disappearance 
of a massive star. We present such a survey, using the first systematic analysis of 
archival Hubble Space Telescope images of nearby galaxies with the aim of identifying 
evolved massive stars which have disappeared, without an accompanying optically 
bright supernova. We consider a sample of 15 galaxies, with at least three epochs of 
Hubble Space Telescope imaging taken between 1994 and 2013. Within this data, we 
find one candidate which is consistent with a 25-30 M© yellow supergiant which has 
undergone an optically dark core-collapse. 
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1 INTRODUCTION 

Massive (>8 M 0 ) stars end their lives as core-collapse su¬ 
pernovae (CCSNe). Once they have evolved off the main 
sequence and developed a Chandrasekhar-mass Fe core, the 
pressure in their core is no longer sufficient to support the 
star against its own gravity. The inner core begins to infall, a 
process which is halted once it reaches nuclear densities and 
forms a proto-neutron star (PNS). At this point the core 
rebounds, and drives a shock outwards through the star. 
According to current simulations, this shock lacks the en¬ 
ergy to halt the stars collapse and explode the outer layers 
of the star. It is believed that the deposition of additional 
energy behind the shock (perhaps in the form of neutrinos) 
can revive it and cause the star to explode as a CCSNe (e.g. 
Burrows 2013; Janka et al. 2007). It is possible, however, for 
a supernova (SN) explosion to fail. If the energy deposition 
is insufficient to power the shock to overcome the infalling 
material, then the shock is not revived (e.g. Woosley 1993; 
Fryer 1999; O’Connor V Ott 2011). This is a more likely 
outcome for more massive stars, which have bigger Fe cores 
and are consequently harder to explode. 

SNe are distinguished primarily by the elements seen in 
their spectra - Type I SNe are H poor, while Type II SNe 
are H rich (Filippenko 1997; Smartt 2009). The Type I 
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SNe are further subdivided into Type lb and Type Ic SNe 
which show the presence or absence of He respectively. H rich 
Type II SNe are separated according to their lightcurves into 
Type IIP (Plateau) SNe, which show an extended period of 
roughly constant luminosity after explosion, and Type IIL 
(Linear) SNe which have a steady decline in luminosity from 
peak. These differing SN types have long been believed to 
result from the extent of the H and/or He envelope of the 
progenitor star at the moment of collapse (Falk V Arnett 
1977), although recently there has been some debate in the 
literature as to whether Type IIP and IIL SNe are distinct 
classes (e.g. Anderson et al. 2014; Faran et al. 2014, and 
references therein). 

Nearby CCSNe present a unique opportunity to test 
this hypothesis by directly identifying their massive stellar 
progenitors in archival images. This was first accomplished 
for SN 1987A in the Large Magellanic Cloud (West et al. 
1987), and for SN 1993J in M51 (Aldering, Humphreys V 
Richmond 1994; Maund V Smartt 2009). Since then, there 
has been significant success in detecting the H rich progeni¬ 
tors of Type II SNe (e.g. Li et al. 2007; Mattila et al. 2008; 
Maund et al. 2011; Van Dyk et al. 2012; Fraser et al. 2014). 
In particular, Type IIP SNe have been shown to come from 
red supergiants (RSGs) with extended atmospheres (Smartt 
et al. 2009). The H-poor progenitors of Type Ibc SNe, how¬ 
ever, have thus far eluded confirmed detections (Eldridge et 
al. 2013), with only iPTF13bvn (Cao et al. 2013; Bersten 


© 0000 RAS 



2 Reynolds, Fraser & Gilmore 


et al. 2014; Eldridge et al. 2015) as a viable, albeit uncon¬ 
firmed, candidate. 

While statistical studies of ensembles of SN progenitors 
(Van Dyk et al. 2003; Smartt et al. 2009) have revealed 
that Type IIP SNe do indeed come from RSGs, there is an 
apparent lack of progenitors with a luminosity comparable 
to that of the brightest known RSGs. Smartt et al. found 
that there was a statistically significant absence of higher 
mass (>16 M 0 ) red supergiants exploding as SNe, terming 
this the “Red Supergiant Problem”. Smartt et al. went on to 
suggest that the red supergiant problem might be explained 
by failed SNe, with stars in the mass range from 16 M 0 to 
30 M 0 collapsing to form black holes with very weak (or 
perhaps non-existent) explosions. 

Furthermore, attempts have been made to constrain SN 
progenitor masses without direct progenitor observations. 
By analysing the resolved stellar population in the vicin¬ 
ity of nearby SNe, the local stellar age can be determined, 
and used to set an upper mass limit on the SN progenitor 
(Williams et al. 2014). Williams et al. set 11 additional (to 
those of (Smartt et al. 2009)) progenitor constraints with 
this method with results consistent with there being no pro¬ 
genitors with >20 M 0 . 

There is some supporting evidence for the suggestion 
that >16 M 0 stars do not explode as optically bright SNe. 
Spectra of SNe taken ^1 yr after explosion can provide a 
diagnostic of the core mass and composition of the star that 
exploded (Mazzali et al. 2010; Jerkstrand et al. 2012). In 
particular, the ejected mass of O, and the strength of the 
[O i] emission lines increase rapidly for a progenitor with 
a mass of more than ~20 M 0 (Jerkstrand et al. 2014). To 
date, however, there has been no observed Type IIP SNe 
with strong [O i] emission consistent with a progenitor with a 
zero-age main sequence (ZAMS) mass >16 M 0 (Jerkstrand 
et al. 2014). We also note that Brown V Woosley (2013) 
found that Galactic abundances can be reproduced even if 
no stars above 25 M 0 contribute metals via SN explosions, 
and that even an even lower threshold of 18 M 0 could be 
accommodated with some modification of uncertain mass 
loss parameters and reaction rates. 

Attempts have also been made to compare the observed 
SN rate with the star formation rate (SFR) (Horiuchi et al. 
2011; Botticella et al. 2012). Horiuchi et al. claimed that 
SN rates were lower than expected given the observed star 
formation rates in galaxies locally and at intermediate red- 
shifts, and suggested that this may be explained by a sig¬ 
nificant number of optically dark (either intrinsically or due 
to dust obscuration) SNe. However, measurements of both 
the SN rate and SFR suffer from signicant systematic and 
random uncertainties (Mathews et al. 2014), and Botticella 
et al. (2012) find that there is no disagreement between the 
SN rate and SFR within 11 Mpc. Mattila et al. (2012) also 
find the SN rate and SFR to be in agreement in the very 
local (<15 Mpc) volume but further argue that a significant 
number of SNe are missed in optical searches beyond the 
local universe due to obscuration by large amounts of dust 
in their host galaxies, especially in luminous and ultralu- 
minous infrared galaxies. Studies including corrections for 
the SNe lost in such highly dust-enshrouded environments 
have found their rates to be well consistent with the SFR 
up to z~1.3 as expected (Dahlen et al. 2012; Melinder et 
al. 2012). 


The relative rates of different CCSN types can also be 
used to probe their progenitor channel. Smith et al. (2011) 
found that 48 per cent of observed CCSNe were of Type 
IIP, and that if these came from the lower extremum of the 
core-collapse mass range, that RSGs between 8.5 and 13.7 
M 0 could account for these. However, Smith et al. went on 
to argue that the “missing” higher mass RSGs may simply 
continue to evolve into a post-RSG phase, before exploding 
as a different type of progenitor. A final piece of evidence for 
the existence of missing SNe lies in the observed distribution 
of compact remnant masses, where Kochanek (2014) argued 
that the discontinuity between the masses of neutron stars 
and black holes could be explained by failed SNe. 

Archival searches for the progenitors of CCSNe such 
as that discussed by Smartt (2009), rely on surveys which 
are searching for bright optical transients. While Lovegrove 
& Woosley (2013) and Piro (2013) have suggested that a 
failed SN may be accompanied by a faint transient, to date 
no convincing candidate has been found. If a core-collapse 
is optically dark (or very faint), then it will not be detected 
by these surveys, and hence its progenitor will not be iden¬ 
tified (although we note that Yang & Lunardini (2011) have 
raised the possibility of detecting failed SNe through their 
neutrino emission). As an alternative approach, Kochanek 
et al. (2008) proposed a “survey about nothing” to monitor 
a sample of nearby galaxies with sufficient depth and reso¬ 
lution to detect individual massive stars. By searching for 
a massive star which had disappeared, it would be possi¬ 
ble to identify faint or dark failed supernovae, which show 
no bright optical display. Kochanek et al. (2008) suggested 
either using the Hubble Space Telescope ( HST ) or ground- 
based 8-m class telescopes for such a survey. Gerke et al. 
(2014) presented the first data from such a survey, using a 
sample of 27 galaxies, observed over 4 years with the Large 
Binocular Telescope. Gerke et al. found a single candidate 
in their data, which they suggested was consistent with a 
failed SN fraction of ~0.3 of all core-collapses. 

In this paper, we present the first attempt to conduct a 
systematic search for failed supernovae using extant archival 
HST data. In Sect. 2 we discuss the selection criteria for 
our sample of galaxies, and the technique used to search for 
disappearing stars; in Sect. 3 we present the candidate failed 
SNe found, while in Sections 4 and 5 we provide a discussion 
of the results and conclusions respectively. 


2 SAMPLE AND METHODOLOGY 

The observational signature of an optically dark core¬ 
collapse SN is a massive star which is seen to disappear 
without an accompanying bright SN. Our basic approach is 
to take galaxies which have been observed with the Hub¬ 
ble Space Telescope (HST) in the same filter, on at least 
three separate occasions, and search for any luminous point 
source which is present in the first two epochs, but is no 
longer present in the third epoch. In the following section 
we elaborate on our methodology, and describe in detail the 
construction of our sample of galaxies and the selection of 
data used. 
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Table 1 . Field of view and spatial resolution of the HST instru¬ 
ments used in this work. 


Instrument 

Field of view 

Pixel scale 

WFPC2 

150" x 150" 

0 . 10 " 

ACS/WFC 

202" x 202" 

0.05" 

WFC3/UVIS 

162" x 162" 

0.04" 

Target and data selection 



The NASA Extragalactic Database (NED) 1 was queried for 
all objects classed as a galaxy; which were listed in either 
the Messier catalog, the New General Catalog, the Index 
Catalog, or the Upsala General Catalog of Galaxies; and 
which had a recessional velocity of less than 2000 km s -1 . 
From this, we obtained a catalog of 2665 galaxies and their 
coordinates, which served as our input sample. 

We adopt a maximum distance limit of 28 Mpc (cor¬ 
responding to a recessional velocity of <2000 km s -1 for 
Ho = 72 km s -1 Mpc -1 , after correcting for infall on the 
Virgo cluster) for our sample. This is the same limit adopted 
by Smartt et al. (2009) for their volume-limited sample of 
Type IIP SN progenitors; beyond this distance many mas¬ 
sive stars are too faint too detect, and the spatial resolution 
of HST causes problems with crowding and blending of point 
sources. 

We consider only data from the three most recent wide- 
held optical cameras on HST , namely the Advanced Camera 
for Surveys / Wide Field Channel (ACS/WFC), the Wide 
Field and Planetary Camera 2 (WFPC2) and the Wide Field 
Camera 3 / Ultraviolet-Visible Channel (WFC3/UVIS). The 
characteristics of each of these instruments are listed in Ta¬ 
ble 1. Images from the High Resolution Channel (HRC) 
of ACS were not used, as the field of view of the camera 
(29 // x25 // ) is significantly smaller. 

As the targets we are searching for are cool RSGs, the 
red F814 W (~ /-band) filter is most sensitive to them. Fur¬ 
thermore, F 8 I 4 W is among the more commonly used HST 
filters, and so data is available for a significant fraction of our 
sample of galaxies. The F 8 I 4 W filters of WFPC2, ACS and 
WFC3 are also quite similar, as can be seen in Fig. 1, allow¬ 
ing us to combine data from multiple instruments. As initial 
criteria for inclusion in our sample, we required a galaxy 
to have been observed on three separate epochs, with at 
least six months between the second and third epochs; and 
with at least some overlapping area common to all three 
epochs. We searched the Mikulski Archive for Space Tele¬ 
scopes (MAST) 2 for any HST observations within a 5 7 ra¬ 
dius around the cataloged position of each of the galaxies in 
our sample, and which met these criteria. The observations 
for each galaxy were then examined using the visualisation 
tools in the Hubble Legacy Archive 3 to determine whether 
their footprints overlapped. 

It is also necessary to have a minimum period between 
the first and second epochs, to ensure that short period 
variables and short-lived transients such as novae can be 
rejected. For most of our sample, we have at least 50 days 


1 http://ned.ipac.caltech.edu/ 

2 http://archive.stsci.edu/ 

3 http://hla.stsci.edu/ 
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Figure 1. Normalized transmission functions for the F 8 I 4 W fil¬ 
ters of WFPC2, ACS and WFC3 onboard HST 

separating the first and second epochs (as listed in Table 
Al), with only one galaxy (NGC 4639) having a shorter pe¬ 
riod. While ideally one would have several years separating 
the first and second epochs, we are limited in this work by 
the available archival data. 

We found 76 galaxies within our input catalog which 
met our initial cuts, and had some overlap between the three 
epochs. For each of these galaxies, NED was queried for any 
non-kinematic measurements of the distance. Kinematic dis¬ 
tances (as used in our initial selection) are prone to error 
for nearby galaxies where peculiar motions will be signifi¬ 
cant compared to recessional velocities. In our analysis of 
candidates, we use Cepheid and SN la luminosity distances 
to the host galaxy. 

We also adopt a minimum distance limit of 10 Mpc for 
our sample. In the case of galaxies closer than this, the field 
of view of images from HST cover only a small fraction of a 
typical galaxy. Furthermore, in a typical 600s exposure with 
HST-\- ACS, a source with magnitude rai=25 is detected at 
~ lOcr, implying that for galaxies within ~ 10 Mpc, many 
sources fainter than Mfsi 4 w < — 5.2 will be detected, which 
are too faint to be massive red supergiants. To avoid galaxies 
with high levels of internal extinction, only those with an in¬ 
clination angle i of <70° were included. As early type galax¬ 
ies have minimal (if any) ongoing star formation, we do not 
expect to find massive red supergiants within them. Galax¬ 
ies with a morphological type of SO or earlier were hence 
excluded from our sample, as were dwarf galaxies (with an 
absolute magnitude Mb < —18) which will contain fewer 
stars than massive spiral galaxies. After making these cuts, 
we were left with a sample of 15 galaxies, which are listed in 
Table 2, along with their basic properties. The imposed cuts 
(distance, inclination angle, galaxy type and galaxy mass) 
should not bias our results in any obvious fashion, but sim¬ 
ply reduce false positives or increase our chance of finding a 
disappearing massive RSG. 

We have also excluded from our sample galaxies cases 
where the data at one of the epochs consisted of a single, 
short (<120 s) exposure. In shallow images such as these, we 
cannot be certain that we do not detect an extant RSG sim¬ 
ply due to low signal to noise. Only three galaxies were ex¬ 
cluded on this basis (NGC 5806, NGC 6217 and NGC 6951). 
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Table 2. Properties of the final galaxy sample in this work. Two of the galaxies have recessional velocities which 
imply a distance of slightly closer than 10 Mpc. However, these should not bias our results because it will not affect 
the detection of massive stars. 


Identifier 

Type 

i (deg) 

v (km s 1 ) 

Distance 61 (Mpc) 

M b (mag) 

Af814W 

M66 

SABb 

68 

788 

6.5 (0.5) 

-21.3 

0.05 

NGC3368 

SABa 

51 

941 

11.8 (0.8) 

-21.0 

0.04 

NGC1566 

SABb 

48 

1221 

19.4 (1.4) 

-21.3 

0.01 

NGC2841 

Sb 

65 

831 

12.3 (0.9) 

-21.0 

0.02 

NGC4496A 

SBd 

39 

1780 

13.6 (1.0) 

-20.2 

0.04 

NGC4321 

SABb 

23 

1677 

14.1 (1.0) 

-22.1 

0.04 

NGC4414 

Sc 

57 

911 

9.03 (0.6) 

-20.2 

0.03 

NGC3982 

SABb 

30 

1354 

21.8 (1.5) 

-19.9 

0.02 

NGC4038 

SBm 

52 

1559 

23.3 (1.6) 

-21.5 

0.07 

NGC4639 

Sbc 

51 

1088 

13.9 (1.0) 

-19.2 

0.04 

NGC4651 

Sc 

50 

913 

14.0 (1.0) 

-19.7 

0.04 

NGC3370 

Sc 

56 

1357 

23.3 (1.6) 

-19.8 

0.05 

NGC3021 

Sbc 

56 

1672 

26.5 (1.9) 

-19.9 

0.02 

NGC1058 

Sc 

58 

634 

9.0 (0.6) 

-18.7 

0.10 

NGC3344 

Sbc 

19 

698 

6.9 (0.5) 

-19.8 

0.05 


a Distances are derived from the recessional velocity from NED, after correction for local motion. When analysing 
our candidates in Sect. 3 we use Cepheid and Type la SN luminosity distances which are discussed in the candidate 
section. 


We also note that cosmic rays, which are typically identified 
using two consecutive exposures, may pose a problem where 
only a single image is available. While there were no such 
cases in our data, in principle cosmic rays can be reliably 
identified in single images using LACOSMIC (van Dokkum. 
2001 ). 

We note that at this point we have not imposed any 
constraints to restrict our sample to galaxies where a SN 
necessarily would have been observed. Within our distance 
limit, systematic surveys such as the Catalina Real Time 
Transient Survey (Drake et al. 2009), the Palomar Transient 
Survey (Law et al. 2009), the Lick Observatory SN Search 
(Leaman et al. 2011), PanSTARRS-1 (Kaiser et al. 2010) 
and ASAS-SN (Shappee et al. 2014) ensure that most SNe 
are detected (although see Prieto et al. 2012 for a counterex¬ 
ample). However, prior to 1998 most of these surveys were 
not in operation, and it is likely that some optically bright 
SNe within this volume were missed (Smartt et al. 2009). 
Another reason why a SN could be missed is due to extinc¬ 
tion, however in such a case the progenitor would almost 
certainly be obscured also. In order to reduce the possibil¬ 
ity that a massive star in our sample exploded as an opti¬ 
cally bright SN, but passed un-noticed, we could exclude any 
galaxies observed prior to 1998. A further requirement could 
be to only consider galaxies which are above +70° or below 
-70° in declination, as such circumpolar targets are visible 
for almost the entire year. This would reduce the chance 
that a SN exploded in a galaxy just as it went behind the 
Sun, and then was too faint to be detected at the start of the 
subsequent observing season 4 . However, as these constraints 
would reduce our sample size still further we have not im- 


4 If one were to restrict the sample to galaxies at high or low 
equatorial latitudes, then one would need to ensure that ground 
based surveys were also monitoring these regions. 


posed them here, but rather discuss their consequences in 
Section 4. 

2.2 Alignment and subtraction 

The data products used in our survey were the Hubble 
Legacy Archive (HLA) pipeline-drizzled mosaics. These have 
been resampled from the original data, and hence are sub- 
optimal for photometry. However, this should not affect dif¬ 
ference imaging to any significant extent. In one case, the 
observations of NGC 1566 from 1995 Aug 30, no mosaic 
was available in the HLA archive; here we created our own 
mosaic using the iraf 5 STSDAS package. 

For each galaxy, before performing difference imaging, 
we registered all frames to a common pixel scale, defined to 
be that of the image taken at the first epoch. To achieve 
this, approximately 30 point sources which were present in 
all three images were identified. As far as possible, these 
sources were distributed evenly across the entire region of 
overlap. We avoided using foreground stars in our alignment, 
as these can show appreciable proper motion in HST data 
taken over timescales of ~years. The centroids of the point 
sources were measured in iraf, and matched coordinate lists 
were used with iraf GEOMAP to derive a geometric trans¬ 
formation between the first and second, and first and third 
epochs. The “general” transform was used in this process, 
which accounts for shifts, rotation, scaling in x and y, and a 
polynomial term which accounts for any residual distortion 
in the images. After deriving a transformation, the images 
from the second and third epochs were then transformed 

5 IRAF is distributed by the National Optical Astronomy Observa¬ 
tory, which is operated by the Association of Universities for Re¬ 
search in Astronomy (AURA) under cooperative agreement with 
the National Science Foundation. 
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into the pixel coordinate system of the image from the first 
epoch using IRAF GEOTRAN. The images were then checked 
visually to ensure that there was no obvious misalignment. 

The different point spread functions (PSFs) of WFPC2, 
ACS and WFC3, in conjunction with the effects of resam¬ 
pling when aligning images, meant that it was not possi¬ 
ble to directly subtract pairs of images from each other. 
Instead, one of the images was first convolved with a kernel 
so that the PSFs of both images matched, before the images 
were scaled to a common flux level and subtracted from 
each other. This technique of “Optimal Image Subtraction” 
(Alard & Lupton 1998) is commonly used in transient sur¬ 
veys to detect SNe, and has been used with HST images to 
confirm the disappearance of progenitor candidates (Maund 
et al. 2014). We used the “High Order Transformation 
of Point Spread Function and Template Subtrac¬ 
tion” (hotpants) package 6 to perform image subtraction 
between the first and third epochs, and second and third 
epoch images for each galaxy in our sample. An example of 
this process is shown in Fig. 2. 

2.3 Source detection and candidate identification 

To identify sources in the difference images, we used the SEX- 
TRACTOR package (Bertin & Arnouts 1996). For a source to 
be detected in the difference images, we require four adja¬ 
cent pixels to be at least 3 a above the measured background. 
For an optically dark SN which is present in the first two 
epochs, and gone in the third epoch, there should be a posi¬ 
tive source in both subtracted images (i.e. Epoch 1 - Epoch 
3, and Epoch 2 - Epoch 3). We hence crossmatched the list 
of detected sources in each difference image, keeping only 
those which were present in both subtracted images, and 
coincident to <4 pixels. Our detection limit varies depend¬ 
ing on the exposure time of the images, and on the local 
background flux at a given position; however, our candi¬ 
dates have F555W magnitudes of <23 - 26.5 mag, and this 
is indicative of the depth to which we are sensitive. 

The source detection and crossmatch process typically 
yielded a list of 50-75 candidates per galaxy. The vast ma¬ 
jority of these candidates are artefacts resulting from the 
subtraction process, such as residuals surrounding bright 
or saturated stars, or residuals in regions of complex back¬ 
ground structure. Examples of such artefacts are shown in 
Fig. 2. To screen these, we performed a visual inspection of 
all candidates, removing bad subtractions, and only keeping 
those sources where there was no obvious source visible in 
the final epoch. At the end of this process we obtained a 
sample of six candidates, which are discussed in detail in 
Sect. 3. 

2.4 Candidate photometry 

For each candidate, we searched the MAST archive for all 
WFPC2, ACS and WFC3 images covering the location of 
the candidate. These images were used in preference to the 
Hubble Legacy Archive products which were used to perform 
difference imaging and identify candidates, as the flux in 

6 http://www.astro.washington.edu/users/becker/v2.0/hotpants 
.html 


these images has not been resampled. The reduced cOf files 
(for WFPC2) and flt or flc files (for ACS and WFC3) 
were downloaded from the MAST archive. 

For the WFPC2 data, PSF-fitting photometry was 
performed with the HSTPHOT package (Dolphin 2000a). 
hstphot has been specifically designed to perform opti¬ 
mally on undersampled WFPC2 images, and includes ze- 
ropoints, charge transfer efficiency (CTE) corrections and 
aperture corrections as per Dolphin (2000b). Before per¬ 
forming photometry, the WFPC2 images were pre-processed 
using the ancillary routines distributed with hstphot. This 
pre-processing included masking pixels which were flagged 
as bad in the data quality frame, combining multiple images 
taken at the same position to reject cosmic rays, masking hot 
pixels, and finally measuring the sky background for each 
pixel by taking the mean of the surrounding pixels. PSF fit¬ 
ting photometry was then performed on the processed im¬ 
ages using pre-computed PSFs appropriate to each filter. 
The most up-to-date CTE corrections were applied to the 
measured photometry for each candidate, along with aper¬ 
ture corrections which were measured from bright stars in 
the images. The magnitudes for each candidate are listed in 
the VEGAMAG system 

For the ACS and WFC3 images, we used the dolphot 
package to perform photometry, dolphot is a modified ver¬ 
sion of hstphot, and includes libraries of PSFs, zeropoints 
and CTE corrections for each filter and instrument. Both 
the ACS and WFC3 images were pre-processed by masking 
bad pixels using the data quality frame for each image, and 
multiplying the image by a pixel area map to account for 
the non-uniform size of pixels due to geometric distortions 
across the instrument field of view. The local sky background 
was then measured for each pixel in the image, based on the 
surrounding pixels. Library PSFs were fitted to sources in 
each frame, and their magnitudes measured. For some of the 
ACS and WFC3 data, flc files were available in the MAST 
archive; these have been corrected for CTE losses in the 
images and so no CTE correction was applied to measured 
magnitudes. For FLT files, CTE corrections were applied to 
the magnitude of each measure source. Aperture corrections 
were applied in all cases. 


3 CANDIDATES 

In the following section, we present the analysis of each of 
the candidates identified. 


3.1 NGC4639-CANDIDATE-1 

The distance to NGC 4639 has been determined using 
Cepheids, mostly recently by Riess et al. (2011) who de¬ 
rived a distance modulus of p=31.67±0.08. This is consis¬ 
tent with the value of p=31.71±0.08 found by Freedman et 
al. (2001), although slightly closer than that obtained from 
the Type la SN 1990N (p = 32.14 ± 0.14; Reindl et al. 
2005). We adopt the Riess et al. distance in the following. 
The foreground extinction towards NGC 4639 is relatively 
low, Ay = 0.071 mag (Schlafly & Finkbeiner 2011). 

The coordinates of NGC4639-CANDIDATE-1 were 
found to be 12 h 42 m 52 s .17 +13°14 / 54 // .8, as measured from 
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Figure 2. Left : The field of view of each of the HST images for NGC 4038/4039 (the Antennae Galaxies), overlaid on an R-band image 
from the CTIO 0.9m telescope (via NED). Right : The results of the image subtraction process for the 2004 and 2010 images of NGC 
4038/4039. The upper panel shows the difference image for the entire galaxy pair, while the lower panels show a zoom-in on (from left 
to right) a bad subtraction from a saturated star, a bad subtraction of a dense knot in a spiral arm, and a good point source in the 
difference image. 
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Figure 3. HST F81fW image cutouts centred on the position of NGC4639-CANDIDATE-1, the location of which is indicated with tick 
marks. All three panels cover a region 100x100 pixels in size, with a common orientation and scale as indicated. The panels show from 
left to right the first, second and third epoch images used to search for candidates; the second and third epoch images have been aligned 
and resampled to match the first epoch as described in the text. 


© 0000 RAS, MNRAS 000 , 000-000 



















Gone without a bang 7 



t y a 


a A 


F300W 

F350LP l—•- 
F439W I—•- 
F547M 
F555W 


V A 


j Ll. 


H - 
H - 


H - 


F814W I # I j 
F160W I • I : 


49800 49900 54100 

MJD 


y a 


I 


55050 


22.5 


CD 

03 


24.5 



23.5 - 


24 - 


53700 


54050 

MJD 


55350 


563f 


Figure 4. The lightcurve of NGC4639-CANDIDATE-1. Circles 
denote WFPC2 observations, triangles are ACS observations and 
squares are WFC3 observations. 


Figure 6. The lightcurve of NGC3021-CAND-1. Triangles denote 
ACS observations, squares are WFC3 observations. Minor tick 
marks are 10 d apart 


Table 3. Photometry for NGC4639-CANDIDATE-1; associated 
errors are given in parentheses. 


Date 

Instrument 

Filter 

Mag 

1995-01-14/15 

WFPC2 

F300W 

22.207 

(0.145) 

- 

- 

F439W 

23.539 

(0.115) 

- 

- 

F547M 

23.446 

(0.093) 

1995-04-24 

- 

F555W 

23.289 

(0.065) 

1995-05-03 

- 

- 

23.248 

(0.052) 

1995-05-10 

- 

- 

23.293 

(0.057) 

1995-05-17 

- 

- 

23.242 

(0.057) 

- 

- 

F814W 

22.814 

(0.041) 

1995-05-21 

- 

F555W 

23.278 

(0.068) 

1995-05-26 

- 

- 

23.321 

(0.076) 

1995-05-30 

- 

- 

23.375 

(0.076) 

1995-06-03 

- 

- 

23.369 

(0.086) 

- 

- 

F814W 

23.003 

(0.064) 

1995-06-06 

- 

F555W 

23.309 

(0.069) 

1995-06-15 

- 

- 

23.406 

(0.071) 

1995-06-23 

- 

- 

23.367 

(0.063) 

1995-07-03 

- 

- 

23.359 

(0.073) 

2006-12-28 

ACS/ WFC1 

F555W 

24.556 

(0.027) 

2007-01-06 

- 

- 

24.520 

(0.025) 

2009-08-07 

WFC3/UVIS1 

F555W 

24.812 

(0.042) 

- 

- 

F814W 

24.671 

(0.091) 

- 

- 

F350LP 

24.780 

(0.038) 


the HLA F814W mosaic. The lightcurve for NGC4639- 
CANDIDATE-1 is plotted in Fig. 4, while selected epochs of 
imaging are shown in Fig. 3. The measured photometry for 
NGC4639-CANDIDATE-1 is listed in Table 3. The F 8 I 4 W- 
band lightcurve shows a 1.5 mag drop between 1995 May 
/ June, and 2009 August, when NGC4639-CANDIDATE- 
1 was detected by dolphot at a much fainter magnitude. 
However, visual inspection of the 2009 images does not re¬ 
veal an obvious point source, and so we cannot exclude the 
possibility that NGC4639-CANDIDATE-1 has disappeared 
completely in F 8 I 4 W in the last epoch, and that dolphot 
is simply detecting some residual background structure at 
its location. 


We calculate the average F 8 I 4 W and F555W magni¬ 
tudes of NGC4639-CANDIDATE-1 over the period April - 
May 1995 to be 23.32 ± 0.05 and 22.91 ± 0.13 mag respec¬ 
tively. Using the adopted extinction and distance towards 
NGC4639, we derive absolute magnitudes of Mfsssw = 
—8.42 d= 0.09 and Mfsi 4 w = —8.80 db 0.15, and a colour 
F555W—F814 W = 0.38 ±0.14 for NGC4639-CANDIDATE- 
1. 

The earliest set of images covering the site of NGC4639- 
CANDIDATE-1 were taken in 1995 using WFPC2, with 
the F218W, F300W , F439W and F547M filters. No source 
was detected in the F218W- filter image. In 1995, NGC4639- 
CANDIDATE-1 was brightest in F300W , suggesting that 
it was relatively hot, and perhaps in outburst at this point. 
Similarly, the robust detection of NGC4639-CANDIDATE-1 
in the F555W filter in 2006 December and 2007 January, at 
approximately 1.2 mags fainter than in 1995, points towards 
a variable source, rather than one which has disappeared en¬ 
tirely by 2009. 


3.2 NGC3021-CANDIDATE-1 

We adopt the Cepheid derived distance modulus of \i = 
32.27 at 0.08 mag for NGC 3021 from Riess et al. (2011). 
The foreground extinction is taken to be My=0.037. 

Using the 2005 HLA F 8 I 4 mosaic, we measured co- 
ordinates of 09 h 50 m 55 s .39 +33°33'14".5 for NGC3021- 
CANDIDATE-1. The region of the candidate is show in Fig. 
5, while the photometry is listed in Table 4 and plotted as 
a lightcurve in Fig. 6. 

From 2005 Nov 12, until 2006 Jan 2, NGC3021- 
CANDIDATE-1 had average magnitudes of F555W=22.98± 
0.03 and F814W—22AA ± 0.04. A year later, over the pe¬ 
riod 2006 November 19-27, the source was only 0.15 mag 
fainter, at F555W = 22.83 ± 0.02 mag. In 2005/2006, 
the source had an absolute magnitudes of —9.33 ± 0.09 
and —9.85 ± 0.09 mag in F555W and F 8 I 4 W respectively, 
given the distance modulus and foreground extinction to¬ 
wards NGC 3021. The F555W-F815W colour of NGC3021- 
CANDIDATE-1 was 0.37 ± 0.04 mag. Some 3.5 years later, 
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Figure 5. HST F814W image cutouts centred on the position of NGC3021-CANDIDATE-1, similar to Fig. 3 


Table 4. Photometry for NGC3021-CANDIDATE-1 


Date 

Instrument 

Filter 

Mag 

2005-11-12 

ACS/WFC2 

F555W 

23.020 

(0.007) 

- 

- 

F814W 

22.500 

(0.006) 

2005-11-22 

- 

F555W 

23.020 

(0.007) 

2005-11-30 

- 

F555W 

22.977 

(0.006) 

- 

- 

F814W 

22.464 

(0.006) 

2005-12-04 

- 

F555W 

22.940 

(0.007) 

2005-12-08 

- 

F555W 

22.936 

(0.006) 

2005-12-11 

- 

F555W 

22.985 

(0.006) 

2005-12-12/13 

- 

F555W 

22.977 

(0.006) 

- 

- 

F814W 

22.442 

(0.006) 

2005-12-14 

- 

F555W 

23.002 

(0.007) 

2005-12-17/18 

- 

F555W 

22.939 

(0.006) 

2005-12-21/22 

- 

F555W 

23.010 

(0.007) 

- 

- 

F814W 

22.407 

(0.006) 

2005-12-26 

- 

F555W 

22.982 

(0.007) 

2006-01-02 

- 

F555W 

22.937 

(0.006) 

- 

- 

F814W 

22.402 

(0.006) 

2006-11-19 

- 

F555W 

22.849 

(0.009) 

2006-11-27 

- 

F555W 

22.819 

(0.009) 

2010-06-03 

WFC3/UVIS2 

F555W 

24.258 

(0.045) 

- 

- 

F814W 

23.503 

(0.059) 

2013-02-07 

WFC3/UVIS1 

F555W 

24.395 

(0.028) 

- 

- 

F814W 

23.690 

(0.041) 


NGC3021-CANDIDATE-1 appears to be ~1.5 mag fainter 
in F555W and ~1 mag fainter in F 8 I 4 W. A final set of 
observations in early 2013 show the source to be similarly 
faint. 

Whether the source detected by dolphot at late times 
is related to NGC3021-CANDIDATE-1 is not immediately 
clear. From examination of the images (Fig. 5), it appears 
that there is some extended background flux at the loca¬ 
tion of the candidate. There may be a contribution to this 
flux from the progenitor, or it may be comprised entirely of 
unrelated background sources. The region of the candidate 
shows several such partially resolved extended sources, so it 
would not be surprising to find a source in close proximity 
to one of these. Alternatively, as massive stars tend to be 
found in star forming complexes, this may be the host clus¬ 
ter of NGC3021-CANDIDATE-1. If this is the case, then 


N 



• 30” 


Figure 7. SDSS r ; image of NGC4496A and NGC4496B. Each 
galaxy is indicated, the dashed ellipses represent the isopho- 
tal major and minor axes at a surface brightness level B =25 
mag/arcsec 2 ). The positions of each of the candidates are marked. 

NGC3021-CANDIDATE-1 clearly dominated the flux of the 
cluster and was its most massive inhabitant, which would 
be consistent with it being the next star to explode. 

3.3 NGC4496-CANDIDATE-1 

NGC4496-CANDIDATE-1 is associated with a galaxy pair, 
NGC 4496A and NGC 4496B, which are also designated 
VV76. While NGC 4496A and NGC 4496B were originally 
thought to be a physically related pair of galaxies, they have 
significantly different recessional velocities of 1700 km s -1 
and 4510 km s -1 respectively (Filippenko et al. 1988). NGC 
4496B, which has a higher recessional velocity has a smaller 
angular size, and is most likely behind NGC 4496A; while 
such coincident superpositions are rare, NGC 4496A and 
NGC 4496B lie toward the Virgo cluster which has a higher 
density of nearby galaxies. 

We measured coordinates of 12 h 31 m 40 s .57 +3°56 / 45 // .l 
for NGC4496-CANDIDATE-1 (Fig. 8). The candidate lies 
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Table 5. Photometry for NGC4496-CANDIDATE-1 


Date 

Instrument 

Filter 

Mag 

1994-05-27 

WFPC2/WF4 

F555W 

26.199 

(0.137) 

1994-06-06 

- 

F555W 

26.485 

(0.266) 

- 

- 

F814W 

23.850 

(0.035) 

1994-06-11 

- 

F555W 

26.251 

(0.134) 

1994-06-16 

- 

F555W 

26.097 

(0.103) 

1994-06-18 

- 

F555W 

26.091 

(0.159) 

1994-06-21 

- 

F555W 

26.132 

(0.149) 

- 

- 

F814W 

23.773 

(0.034) 

1994-06-22 

- 

F555W 

25.966 

(0.160) 

1994-06-26 

- 

F555W 

26.128 

(0.169) 

1994-06-29 

- 

F555W 

26.128 

(0.169) 

1994-07-01 

- 

F555W 

26.152 

(0.157) 

1994-07-03 

- 

F555W 

26.260 

(0.120) 

1994-07-16 

- 

F555W 

26.239 

(0.236) 

- 

- 

F814W 

23.829 

(0.043) 

1994-07-20 

- 

F555W 

26.258 

(0.172) 

1994-07-26 

- 

F555W 

26.014 

(0.167) 

1994-08-02 

- 

F555W 

26.053 

(0.025) 

1994-08-06 

- 

F555W 

26.355 

(0.259) 

- 

- 

F814W 

23.833 

(0.043) 

1994-08-07 

- 

F555W 

26.289 

(0.193) 

2001-03-17 

WFPC2/WF2 

F814W 

>24.9 



49500 52000 

MJD 


Figure 9. The lightcurve of NGC4496-CANDIDATE-1. Circles 
denote WFPC2 observations, minor tick marks are 10 d apart. 

within the isophotal radius of NGC 4496A at B — 25, and 
well outside the isophotal radius of NGC 4496B. NGC4496- 
CANDIDATE-1 also lies to the north of NGC 4496A, 
whereas NGC 4496B lies to the south. It hence appears al¬ 
most certain that NGC4496-CANDIDATE-1 is associated 
with NGC 4496A. The distance towards NGC 4496A has 
been measured both using Cepheids, and with the peak lu¬ 
minosity of the Type la SN 1960F. We adopt a Cepheid- 
based distance modulus of 30.86 d= 0.03 mag, corresponding 
to 14.86 Mpc (from the final distance modulus with period 
cuts and correction for metallicity; Freedman et al. 2001). 
This is consistent with the luminosity distance from SN 
1960F, 30.86 mag (Reindl et al. 2005). We adopt a fore¬ 
ground extinction of Ay = 0.067 mag towards both NGC 
4496A and NGC 4496B. 


The measured magnitudes for NGC4496- 
CANDIDATE-1 are fisted in Table 5, and plotted as 
a lightcurve in Fig. 9. The source appears to have a 
constant magnitude in both F555W and F814 W over the 71 
days it was observed in mid- 1994, with average magnitudes 
of 26.18 ± 0.13 and 23.82 zb 0.03 mag in F555W and F814W 
respectively. The quoted uncertainties on the average 
magnitude are the standard deviation of all measurements, 
which in both cases are comparable to the photometric 
uncertainty on individual measurements. 

The absolute magnitude of NGC4496-CANDIDATE-1 
in 1994 is relatively faint for a massive RSG. After cor¬ 
recting for distance and foreground extinction, the source 
has an absolute magnitude of F555W— —4.75 zb 0.13 mag, 
and a F555W—F814W colour of 2.33 mag. Such a magni¬ 
tude is similar to that of the faintest RSGs, which have 
been claimed as sub-luminous Type IIP progenitors (Fraser 
et al. 2011). Nearly seven years later, in 2001 March, the 
source is no longer detected in F814 W images. From sur¬ 
rounding sources, we set a 5cr limit of F814W> 24.9 mag on 
NGC4496-CANDIDATE-1, implying it has faded by a least 
a magnitude. 

3.4 NGC4496-CANDIDATE-2 

NGC4496-CANDIDATE-2 (Fig. 10) was also found in the 
NGC 4496A/B pair. Similarly to NGC4496-CANDIDATE- 
1, NGC4496-CANDIDATE-2 lies closer to NGC4496A, is 
within its isophotal radius, and is outside the radius of NGC 
4496B. We adopt the same distance modus and foreground 
extinction (|i = 30.86 zb 0.03, and Ay = 0.067 mag) as 
used for NGC4496-CANDIDATE-1. We show the lightcurve 
of NGC4496-CANDIDATE-2 in Fig. 11, and report the 
magnitude of the source in each frame in Table 6. The 
coordinates of NGC4496-CANDIDATE-2 are 12 h 31 m 41 s .07 
+3°56 / 17 // .3. 

NGC4496-CANDIDATE-2 is either not detected, or 
marginally detected in F555W in 1994. Over the first set of 
images, the average magnitude of NGC4496-CANDIDATE- 
2 is F555W = 26.46 ±0.18, F 8 I 4 W = 23.93±0.04 mag. This 
implies an absolute magnitude of F555W = —4.47±0.18 and 
F 8 I 4 W = —6.97 ± 0.05 mag. We note that in many of the 
F555W images the source is not detected; this is not sur¬ 
prising, as the source is quite faint in this filter. The source 
is detected in all F555W images when there is contempora¬ 
neous F 8 I 4 W data, as we can use the latter to determine 
the position of the source in the former (effectively carrying 
out forced photometry at the position of the source in the 
F555W image). As in the case of NGC4496-CANDIDATE- 
1, if NGC4496-CANDIDATE-2 is a RSG, then it would be 
relatively faint, and hence low mass. 

An alternative interpretation for NGC4496- 
CANDIDATE-2 is that it is in fact associated with 
NGC 4496B to the south. The distance to NGC 4496B 
based on its recessional velocity (after correction for infall 
on Virgo) is 70.1 Mpc, which gives a distance modulus 
of 34.23 mag. This would make the absolute magnitude 
of NGC4496-CANDIDATE-2 -7.84 mag in F555W and 
— 10.34 mag in F 8 I 4 W. However, at this distance the 
implied galacto-centric radius of NGC4496-CANDIDATE-2 
would be 14.6 kpc, which is quite large. On this basis, and 
on the fact that NGC4496-CANDIDATE-2 lies well outside 
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Figure 8. HST F 8 I 4 W image cutouts centred on the position of NGC4496-CANDIDATE-1, similar to Fig. 3. 
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Figure 10. HST F 8 I 4 W image cutouts centred on the position of NGC4496-CANDIDATE-2, similar to Fig. 3. 
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Figure 11. The lightcurve of NGC4496-CANDIDATE-2. Circles 
denote WFPC2 observations, minor tick marks are 10 d apart. 


the isophotal radius of NGC 4496B, we regard it as more 
likely to be associated with NGC 4496A. 


3.5 NGC4496-CANDIDATE-3 

The coordinates of NGC4496-CANDIDATE-3 (Fig. 12) are 
12 h 31 m 40 s .72 +3°55 , 57 ,, .2. The source lies equidistant be¬ 
tween NGC 4496A and NGC 4496B, and is within the B— 25 
isophote of both galaxies, making it difficult to conclusively 
associate it with either. If we take the B band isophote as 
measuring the outer edge of the galaxy, then 65 per cent 
of the area of NGC 4496B is interior to the position of 
NGC4496-CANDIDATE-3, compared to only 16 per cent 
of the area of NGC 4496A. 

The photometry for NGC4496-CANDIDATE-3 is listed 
in Tab. 7, while the lightcurve is shown in Fig. 13. The 
source appears to show variability over the ~70 day pe¬ 
riod of observations in 1994, as it increases in brightness 
in F 8 I 4 W by 0.27 mag. The F555 W lightcurve shows sim¬ 
ilar structure, with an increase over the period of observa¬ 
tions in 1994, with an amplitude of ~0.5 mag. The average 
F555W and F814 mags are 25.45 ± 0.10 and 23.50 ± 0.11 
mag respectively, corresponding to an absolute magnitude 
of F555W — —5.48 =t 0.11 if associated with NGC 4496A, or 
—8.85 mag if it lies in the more distant host NGC 4496B. 
In either case, the F555W—F814W colour of NGC4496- 
CANDIDATE-3 is 1.92 ± 0.15 mag. 

As for the other candidates in NGC 4496A/B, the ab- 
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Table 6. Photometry for NGC4496-CANDIDATE-2 Table 7. Photometry for NGC4496-CANDIDATE-3 


Date 

Instrument 

Filter 

Mag 

Date 

Instrument 

Filter 

Mag 

1994-05-27 

WFPC2/PC 

F555W 

<27.3 

1994-05-27 

WFPC2/PC 

F555W 

25.778 (0.073) 

1994-06-06 

- 

F555W 

26.616 (0.285) 

1994-06-06 

- 

F555W 

25.552 (0.063) 

- 

- 

F814W 

23.938 (0.045) 


- 

F814W 

23.628 (0.038) 

1994-06-11 

- 

F555W 

<27.3 

1994-06-11 

- 

F555W 

25.480 (0.055) 

1994-06-16 

- 

F555W 

<27.3 

1994-06-16 

- 

F555W 

25.435 (0.059) 

1994-06-18 

- 

F555W 

<27.3 

1994-06-18 

- 

F555W 

25.393 (0.060) 

1994-06-21 

- 

F555W 

26.187 (0.200) 

1994-06-21 

- 

F555W 

25.512 (0.091) 

- 

- 

F814W 

23.887 (0.041) 


- 

F814W 

23.525 (0.029) 

1994-06-22 

- 

F555W 

26.293 (0.108) 

1994-06-22 

- 

F555W 

25.483 (0.062) 

1994-06-26 

- 

F555W 

<27.3 

1994-06-26 

- 

F555W 

25.406 (0.074) 

1994-06-29 

- 

F555W 

<27.3 

1994-06-29 

- 

F555W 

25.383 (0.060) 

1994-07-01 

- 

F555W 

26.595 (0.144) 

1994-07-01 

- 

F555W 

25.362 (0.047) 

1994-07-03 

- 

F555W 

<27.3 

1994-07-03 

- 

F555W 

25.335 (0.061) 

1994-07-16 

- 

F555W 

26.598 (0.344) 

1994-07-16 

- 

F555W 

25.477 (0.065) 

- 

- 

F814W 

23.983 (0.051) 


- 

F814W 

23.486 (0.040) 

1994-07-20 

- 

F555W 

<27.3 

1994-07-20 

- 

F555W 

25.321 (0.064) 

1994-07-26 

- 

F555W 

<27.3 

1994-07-26 

- 

F555W 

25.443 (0.068) 

1994-08-02 

- 

F555W 

<27.3 

1994-08-02 

- 

F555W 

25.433 (0.069) 

1994-08-06 

- 

F555W 

26.460 (0.241) 

1994-08-06 

- 

F555W 

25.406 (0.067) 

- 

- 

F814W 

23.930 (0.053) 


- 

F814W 

23.356 (0.027) 

1994-08-07 

- 

F555W 

<27.3 

1994-08-07 

- 

F555W 

25.452 (0.065) 

2001-03-17 

WFPC2/WF4 

F814W 

<25.000 

2001-03-17 

WFPC2/WF4 

F814W 

<25 


solute magnitude of NGC4496-CANDIDATE-3 is relatively 
faint for a massive RSG if it lies in NGC 4496A. However, 
in this case the source could be plausibly associated with 
NGC 4496B. If associated with NGC 4496B, then NGC4496- 
CANDIDATE-3 has an absolute magnitude of F555W~ 
—9. On a Hertzsprung-Russell (HR) diagram, NGC4496- 
CANDIDATE-3 lies between the evolutionary tracks of a 
25 and a 30 M© star, as shown in Fig. 16. The colour of 
NGC4496-CANDIDATE-3 is also consistent with that ex¬ 
pected for an M-type supergiant. On this basis, NGC4496- 
CANDIDATE-3 is consistent with a massive RSG, if associ¬ 
ated with NGC 4496B. On the basis of the available data, it 
is impossible to determine whether this is indeed the case. 


3.6 NGC4321-CANDID ATE-1 

NGC4321-C ANDIDATE-1 is located at 12 h 22 m 58 3 .37 
+15°47'58".0 (Fig. 14) in NGC 4321 (also known as M100). 
We adopt a distance modulus for NGC 4321 of p = 30.91 ± 
0.07 (Freedman et al. 2001), and a foreground extinction of 
Ay = 0.072 mag. The measured magnitudes for NGC4321- 
CANDIDATE-1 are listed in Table 8, while a lightcurve is 
plotted in Fig. 15. 

The source was first observed over 57 days from April 
to June in 1994, when it had an average F814W magni¬ 
tude of 23.40 mag. During this period the source appears 
to be variable, with an amplitude of 0.26 mag in F 8 I 4 W , 
which is considerably larger than the typical photometric 
error of ~0.03 mag. The variability in F555W is even larger 
at 0.76 mag, although the F555W photometry has larger 
uncertainties as the source is fainter in this band. Two years 
later in 1996 April, the source had a similar magnitude in 
both F555W and F 8 I 4 W. However, by 1994 July, the source 
had faded by ~0.7 mag in the F 8 I 4 W band. A final set of 
observations were taken much later in 2008 January, and 



Figure 13. The lightcurve of NGC4496-CANDIDATE-3. Circles 
denote WFPC2 observations, minor tick marks are 10 d apart. 

while the F 8 I 4 W filter was not used, an F791W filter im¬ 
age was obtained. F791W has a fairly similar bandpass to 
F 8 I 4 W 1 apart from having a red cut-off which is ~500A 
shorter than that of the latter. We measured a F791W mag¬ 
nitude of 24.058 zb 0.158, which is identical to the F 8 I 4 W 
magnitude from 1996 July. 

The absolute magnitude of NGC4321-CANDIDATE-1 
is F555W = —5.26, and the F555W—F814 colour is 2.29 db 
0.22 during 1994. 

Based on its photometric evolution, it appears that 
NGC4321-CANDIDATE-1 is a variable source, rather 
than a star which has undergone a terminal event. 
From the relatively sparsely sampled lightcurve, NGC4321- 
CANDIDATE-1 seems to have transitioned from a “bright” 
state in 1994, to a “faint” state over the following years. 
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Figure 12. ThST F814W image cutouts centred on the position of NGC4496-CANDIDATE-3, similar to Fig. 3. 
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Figure 14. HST F 8 I 4 W image cutouts centred on the position of NGC4321-CANDIDATE-1, similar to Fig. 3. 
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Figure 15. The lightcurve of NGC4321-CANDIDATE-1. Circles 
denote WFPC2 observations, minor tick marks are 10 d apart 
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3.7 Candidates summary 

Of the six candidates for which we conducted further 
analysis, both NGC4639-CANDIDATE-1 and NGC4321- 


CANDIDATE-1 appeared to be variable. NGC4496- 
CANDIDATE-1 and NGC4496-CANDIDATE-2 were too 
faint to be consistent with massive RSGs, while NGC4496- 
CANDIDATE-3 is only a plausible massive star if associated 
with the more distant galaxy in the NGC 4496A/B pair. We 
hence find only a single source, NGC3021-CANDIDATE-1, 
which remains as a good candidate for a massive star which 
underwent a failed SN. 


4 DISCUSSION 
4.1 Expected numbers 

To determine rates (or a limit to the rate) of failed SNe from 
the data presented in this paper is a non-trivial task. Ideally, 
one would quantify the detection efficiency for disappearing 
sources at all positions in each galaxy in our sample using 
artificial star tests. The results of this will vary depending 
on the magnitude of the source, the quality of the data, and 
also on the local background and degree of crowding. Once 
detection efficiencies have been computed, one can weight 
the detection probability at each position in the galaxy by 
the local star formation rate, and the initial mass function 
(and hence progenitor magnitude) given the expected mass 
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Table 8. Photometry for NGC4321-CANDIDATE-1 


Date 

Instrument 

Filter 

Mag 

1994-04-23 

WFPC2/WF2 

F555W 

25.500 

(0.137) 

- 

- 

F814W 

23.258 

(0.031) 

1994-05-04 

- 

F555W 

25.203 

(0.160) 

1994-05-06 

- 

F555W 

25.682 

(0.140) 

1994-05-09 

- 

F555W 

25.854 

(0.017) 

1994-05-12 

- 

F555W 

25.793 

(0.232) 

- 

- 

F814W 

23.386 

(0.029) 

1994-05-16 

- 

F555W 

25.758 

(0.138) 

1994-05-20 

- 

F555W 

25.863 

(0.103) 

1994-05-26 

- 

F555W 

25.781 

(0.113) 

1994-05-31 

- 

F555W 

25.753 

(0.133) 

1994-06-07 

- 

F555W 

25.743 

(0.180) 

- 

- 

F814W 

23.422 

(0.043) 

1994-06-17 

- 

F555W 

25.728 

(0.017) 

- 

- 

F555W 

25.960 

(0.325) 

1994-06-19 

- 

F814W 

23.521 

(0.047) 

- 

- 

F555W 

25.769 

(0.201) 

1996-04-27 

- 

F814W 

23.343 

(0.037) 

1996-07-09 

WFPC2/PC 

F336W 



- 

- 

F439W 



- 

- 

F555W 



- 

- 

F675W 

25.419 (0.208) 

- 

- 

F814W 

24.058 (0.096) 

2008-01-04 

WFPC2/WF3 

F380W 



- 

- 

F439W 



- 

- 

F555W 

25.885 (0.170) 

- 

- 

F702W 

24.723 (0.128) 

- 

- 

F791W 

24.058 (0.158) 


range of failed SNe. However, such a calculation would ne¬ 
cessitate numerous assumptions (such as the unknown mag¬ 
nitude range of failed SN progenitors). 

As a simpler alternative, which provides a useful sanity 
check on our results, we simply ask the question “Is the 
number of failed SN candidates seen in our sample consistent 
with the number of optically bright CCSNe seen in the same 
volume and timespan?” This approach naturally accounts 
for unknowns such as the local star formation rate at any 
given location, and the irregular area of overlap between 
multiple HST fields of view. We find that of our sample, 
only SNe 2009hd and 2004gt lie within the area common to 
all our images, and exploded between the second and third 
epoch of imaging. If these two SNe had come from bright 
RSGs, which disappeared completely, then we should have 
detected them in our sample. However, neither SN 2004gt or 
SN 2009hd were detected as failed candidates in our data. 
SN 2004gt was a Type Ibc SN, and as deep searches for 
its progenitor in the same archival HST images as we have 
used were unsuccessful (Gal-Yam et al. 2005; Maund et al. 
2005), it is unsurprising that we do not detect it as a failed 
candidate. For SN 2009hd (Elias-Rosa et al. (2011)), the 
third epoch images we used were taken of SN 2009hd itself, 
and as the SN is still bright in these we cannot test for the 
disappearance of its progenitor. 

If we assume that approximately ~25 per cent of mas¬ 
sive stars may end their lives as optically dark SNe (based 
on the mass range for Type II SN progenitors inferred by 
Smartt et al. 2009), and two bright SNe fell within our sam¬ 
ple, then we may expect to see around ~0.7 failed SNe 


in our sample. As discussed in Sect. 3, only NGC3021- 
CANDIDATE-1 remains as a plausible failed SN candidate 
(along with NGC4496-CANDIDATE-3 if we adopt the larger 
distance for it). We hence regard the number of candidate 
failed SNe found as being broadly consistent with expecta¬ 
tion, albeit with large uncertainties given the small sample 
size. If instead we use the results of Gerke et al. (2014), 
who found that 30^23 P er cent of core-collapses resulted in 
failed SNe, we expect to find around ~0.9 failed SNe in our 
sample. This is again broadly consistent, albeit with uncer¬ 
tainties too large to say anything substantive about failed 
SN rates. 


4.2 What do optically dark SNe really look like? 

While the premise of this paper was to search for the dis¬ 
appearance of a massive and luminous RSG without an ac¬ 
companying optical display, it is worth considering in more 
detail whether any transient would be associated with such 
an event. There have been some attempts to answer this 
question from a theoretical basis (e.g. Nadezhin 1980; Love- 
grove V Woosley 2013; Piro 2013), which have focused on 
the effect of removing around O.IM© of mass from the core of 
a massive star via neutrino emission. The outer layers of the 
star will adjust to this rapid loss of mass, and in the models 
of Lovegrove V Woosley (2013) a fraction of the envelope is 
ejected with a kinetic energy of about 10 47 ergs. This would 
be accompanied by a faint red transient which lasts some 
hundreds of days. Piro (2013) suggested that there would 
be shock breakout emission from such an event, peaking in 
the ultraviolet (UV) and lasting only a few days. 

A transient with brief, UV bright emission would likely 
not be seen in our data, as our typical cadence is extremely 
low. However, a slow red transient would most likely be 
seen in our data. Lovegrove V Woosley (2013) calculated 
that such a transient would have a luminosity of 10 39 - 
10 4 ° ergs s -1 . At a typical distance of ^20 Mpc, this would 
have an apparent magnitude of 20. While such a tran¬ 
sient may well be missed in ground based surveys with 1-m 
or 2-m class telescopes, which typically have a limiting mag¬ 
nitude of > 21 (or even shallower in the case of the LOSS 
survey, which had a limiting magnitude of 18-19.5 mag; Lea- 
man et al. 2011), it would be clearly visible in our HST data. 

4.3 Contaminants 

A major concern for any attempt to find optically dark SNe 
is that variable stars or other transients could mimic the 
disappearance of a massive star, either through appearing 
brighter in the first two epochs of imaging, or fainter in the 
third epoch. 

R Coronae Borealis (RCB) stars are one such possi¬ 
ble contaminant. These are evolved, hydrogen poor F or G 
supergiants, which experience repeating, irregular episodes 
where they decline by up to 8 magnitudes for hundreds of 
days (Clayton 1996). This is thought to be caused by the 
formation of clouds of amorphous carbon dust which can 
temporarily obscure the star. Related to RCBs are DY Per 
type stars, which are cooler than typical RCBs with effective 
temperatures of ~3500 K. RCBs are quite rare, with only 
~50 known in the Milky Way, while DY Per type stars are 
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Figure 16. A colour magnitude diagram for the candidates presented in this paper (black points), compared to solar metallicity stellar 
evolutionary tracks from STARS models (Eldridge & Stanway 2009). The plotted magnitude for each source is the average magnitude over 
the first two epochs of imaging, corrected for distance and foreground reddening. The zero-age main sequence mass of each evolutionary 
model is indicated. The erratic behaviour of the tracks in the lower left is due to the more massive stars entering the Wolf-Rayet phase. 
The F555W and F814W magnitudes for the candidates are plotted, while the evolutionary models are in the Johnson-Cousins system, 
however these are sufficiently similar (with differences of a few hundredths of a magnitude) to permit a direct comparison.Two possible 
positions are shown for NGC4496-CANDIDATE-1, depending on whether it is associated with NGC 4496A or NGC 4496B. The green 
shaded region shows the range of RCB star magnitudes from Tisserand et al. (2009). 


even less common. The absolute magnitude of RCBs are, 
however, fainter than the candidates found in our sample. 
The region occupied by known RCB stars from Tisserand et 
al. (2009) is shown in Fig. 16, and it is clear that all of our 
candidates lie above this region. We hence discount RCBs 
as a possible explanation for our candidates. 

Pulsating stars, which lie in the instability strip of the 
HR diagram are another source of concern. Some, such as 
the RV Tauri variables, can be dismissed as being too faint 
to appear in our data; even the brightest R Scuti has a lu¬ 
minosity log L/L 0 <4 dex. However, some long period vari¬ 
ables, and in particular Mira variables, are more of a con¬ 
cern. Miras are pulsating, late type stars which can have 
amplitudes of several magnitudes, and periods of hundreds 
of days. To assess whether any of our candidates could be 
Miras, in Fig. 17 we compare them to a sample of Miras 
found in the Large Magellenic Cloud by the OGLE survey 
(Soszynski et al. 2009; Udalski et al. 2008). The population 
of Mira type variables in the LMC appears to be fainter, and 
for the most part significantly redder than our candidates. 

Finally, we consider the possibility that known classes 
of transients could mimic a disappearing RSG. This would 



V-l 


Figure 17. Colour magnitude diagram for Mira variables dis¬ 
covered in the LMC by the OGLE survey (Soszynski et al. 2009; 
Udalski et al. 2008). The mean /-band magnitudes of the OGLE 
Miras are plotted as grey points, with error bars corresponding to 
their amplitude. The failed SN candidates discussed in this paper 
are shown as black diamonds. 
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be more of a concern if we had only two detections before 
a source is seen to disappear. In such a scenario, the first 
and second detections of a source could be two successive 
outbursts, perhaps of a recurrent nova, or of an S-Doradus 
type luminous blue variable, which happened to be observed 
at a similar magnitude. In the third image, the hypotheti¬ 
cal source could appear to have disappeared as it is not in 
outburst. However, as each of the candidates in our sam¬ 
ple are detected on more than six separate occasions before 
disappearing, we regard this as unlikely. 


4.4 Was an optically dark SN seen in our sample? 

We cannot conclusively say whether an optically dark SN 
was seen in our sample. Our best candidate, NGC3021- 
CANDIDATE-1 appears to lie in a complex of bright sources 
in a star forming arm of a massive spiral galaxy. On the HR 
diagram shown in Fig. 16, NGC3021-CANDIDATE-1 lies 
close to the track for a 30M© star. The candidate lies above 
the 25M© track, and so we adopt a mass of 25-30 M 0 . In 
the STARS models plotted, if the star lies at the upper end 
of this mass range, then it may not end its life as an RSG, 
but rather cross the HR diagram to the blue to become a 
Wolf Rayet star. 

The colour of NGC3021-CANDIDATE-1 is bluer than 
that of an RSG, and is in fact consistent with an F8 super¬ 
giant (Pickles 1998). It is important to note however, that 
the effective temperature of a star at the point it explodes 
is dependent on its radius and the extent of any mass loss. 
In the case of SN 201 ldh, a yellow supergiant (YSG) was 
observed to explode as a core-collapse SN (Maund et al. 
2011; Van Dyk et al. 2011), and while the progenitor of SN 
201 ldh likely had its envelope stripped by a binary compan¬ 
ion (Folatelli et al. 2014), it at least establishes a precedent 
for massive stars exploding as YSGs. Regardless of the pro¬ 
genitor temperature, the progenitor luminosity should scale 
monotonically with core (and hence progenitor) mass. 

We checked all available online databases for any 
reported SN or transient at the position of NGC3021- 
CANDIDATE-1 between 2006 and 2013. NGC 3021 was tar¬ 
geted by the LOSS survey (Leaman et al. 2011), and was also 
within the footprint of PanSTARRS-1, PTF and CRTS. To 
the best of our knowledge, none of these surveys reported 
a SN or a candidate transient at the location of NGC3021- 
CANDID ATE-1. 

While NGC4496-CANDIDATE-3 is consistent with a 
massive RSG if associated with the more distant galaxy 
NGC 4496B, it is impossible to determine this on the basis 
of the available data. Furthermore, if associated with NGC 
4496B, then the greater distance of the host would formally 
exclude this from our volume limited sample. 


5 CONCLUSIONS 

HST imaging has been shown to be a powerful tool in the 
search for disappearing massive stars. In our preliminary 
study we have found one candidate failed SNe which mer¬ 
its further investigation. HST delivers an order of magni¬ 
tude improvement in spatial resolution over ground based 
data. We note that many of our candidates would have been 


blended with nearby sources in ground based images, ren¬ 
dering their identification more difficult. 

With the successful demonstration of image subtraction 
using archival HST images, one can begin to consider design¬ 
ing an experiment to increase the number of candidates seen. 
We note that while we required galaxies which had been ob¬ 
served three times with HST , there are over one hundred 
more galaxies which had two available epochs of imaging. 
Approximately 30 of these have some overlap between the 
first and second epoch, hence if one were to obtain a third 
epoch of tailored observations for each of these, our sample 
size would be increased by a factor of >2. 

We estimate that to monitor a galaxy would require 
an image every 4 months, ideally in two filters ( F555W and 
F 8 I 4 W). A 4 month cadence would ensure that faint, but rel¬ 
atively long-lived transients such as those modelled by Love- 
grove k Woosley (2013) would not be missed, and would also 
permit long period variable stars such as Miras to be identi¬ 
fied. In a single orbit with HST -\-ACS one can reach limiting 
magnitudes of V = 27.9 and I = 27.4; at a distance of 20 
Mpc this corresponds to a limiting absolute magnitude of 
V — —3.6 and I — —4.1. So, to monitor a sample of 20 
galaxies to this depth one would require 120 orbits of HST 
time per year, a significant investment of time. The prob¬ 
ability of success could be increased by selecting galaxies 
which have a relatively prolific SN rate, such as NGC 4303 
which has hosted four CCSNe in the last 16 years. Such a 
survey would also facilitate a wide range of science, includ¬ 
ing monitoring the variability of massive stars, and searching 
for a range of other faint transients such as stellar mergers, 
common-envelope ejections and luminous red novae (i.e. the 
“gap transients” which lie between novae and supernovae; 
2007 Kulkarni k Kasliwal). 

Our original survey was designed to find failed SNe with 
only three epochs of data - two early epochs in which a star 
was present, and a late epoch in which the same star had 
vanished. However, the candidates presented in this work 
appear to have been preferentially identified in galaxies for 
which significantly more data were available. This suggests 
that the original expectation that three epochs would be suf¬ 
ficient to find failed SNe may be overly optimistic, although 
with deep followup observations we could still verify many 
candidates. Running an ongoing HST survey with higher 
cadence monitoring would help to separate out failed SNe 
from contaminants. 

Finally, we note that alongside the ongoing observa¬ 
tional project to find failed SNe and their progenitors, there 
is a pressing need for theoretical work on the mass ranges 
of stars which are expected to explode, and to understand 
where the islands of “explodability” (e.g. Ugliano et al. 
(2012); Pejcha k Thompson (2014); Ertl et al. (2015)) may 
lie. 
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Table Al. Archival images used for each galaxy. At is the period (i 
of that galaxy. 



Date 

Instrument 

Exposure (s) 

At (yr 

NGC 1058 

Epoch 1: 

2001-07-03 

WFPC2 

460 

- 

Epoch 2: 

2007-10-29 

WFPC2 

320 

6.33 

Epoch 3: 

2008-11-20 

WFPC2 

700 

7.39 

NGC 1566 

Epoch 1: 

1995-08-30 

WFPC2 

300 

- 

Epoch 2: 

2009-03-11 

WFPC2 

600 

13.54 

Epoch 3: 

2013-09-03 

WFC3 

989 

18.02 

NGC 2841 
Epoch 1: 

2000-02-29 

— 2000-04-09 

WFPC2 

44000 


Epoch 2: 

2005-04-05 

WFPC2 

500 

4.6 

Epoch 3: 

2010-01-05 

WFC3 

1412 

9.36 

NGC 3021 

Epoch 1: 

2005-11-12 

ACS/WFC 

4800 

- 

Epoch 2: 

2006-01-02 

ACS/WFC 

4800 

0.14 

Epoch 3: 

2013-02-07 

WFC3 

1126 

7.64 

NGC 3344 

Epoch 1: 

2001-11-13 

WFPC2 

460 

- 

Epoch 2: 

2012-02-24 

ACS/WFC 

900 

10.29 

Epoch 3: 

2014-03-02 

WFC3 

980 

12.31 

NGC 3368 - 
Epoch 1: 

Position 

1994-04-06 

WFPC2 

3000 


Epoch 2: 

1995-04-06 

WFPC2 

3000 

1.00 

Epoch 4: 

2009-11-12 

WFC3 

350 

15.60 

NGC 3368 - 
Epoch 1: 

Position 2 

2000-12-12 

WFPC2 

700 


Epoch 2: 

2001-11-06 

WFPC2 

320 

0.90 

Epoch 4: 

2009-11-12 

WFC3 

350 

8.93 

NGC 3370 

Epoch 1: 

2001-07-06 

WFPC2 

460 

- 

Epoch 2: 

2003-03-31 

- 2003-05-22 

ACS/WFC 

24000 

1.73 

Epoch 3: 

2013-01-06 

WFC3 

1136 

11.50 

NGC 3627 

Epoch 1: 

2001-02-24 

WFPC2 

700 

- 

Epoch 2: 

2001-05-28 

WFPC2 

700 

0.25 

Epoch 3: 

2009-12-14 

ACS 

520 

8.81 

NGC 3982 

Epoch 1: 

2000-04-02 

WFPC2 

5000 

- 

Epoch 2: 

2001-10-06 

WFPC2 

460 

1.51 

Epoch 3: 

2009-12-14 

WFC3 

700 

9.35 

1 Mosaic created manually as 

discussed in 

the text. 



years) between a given epoch and the earliest observations 


) HLAfilename or Dataset 

hst_09042_n2_wfpc2_f814w_wf_drz.fits 
hst _11119_01 _wfpc2 _f814w_wf _drz. fits 
hst_10877_32_wfpc2_f814w_wf_drz.fits 


1995-08-30_mosaic_NGC1566.fits 1 
hst_11987_15_wfpc2_f814w_wf_drz.fits 
hst _13364_17_wfc3 _uvis_f814w_drz .fits 


hlsp_sgal_hst_wfpc2_n2841_f814w_vl_mosaic-sci_drz.fits 

hst_10402_12_wfpc2_f814w_wf_drz.fits 

hst_11360_r2_wfc3_uvis_f814w_drz.fits 


HST_10497_13_ACS_WFC_F814W_drz.fits 
HST_10497_24_ACS_WFC_F814W_drz.fits 
hst _12880_76_wfc3_uvis_f814w_drz.fits 


hst_09042_40_wfpc2_f814w_wf_drz.fits 
j bt 10k020_drz .fits 
hst _13364_22_wfc3_uvis_f814w_drz.fits 


hst_05415_01_wfpc2_f814w_wf.fits 
hst _05415_0p _wfp c 2 _f814 w _wf. fits 
hst _11646_02 _wfc3 _u vis_f814w_drz. fits 


hst _08602_56 _wfp c2 _f814 w _wf _drz .fits 
hst_09042_n5_wfpc2_f814w_wf_drz.fits 
hst _11646_02 _wfc3 _uvis_f814w_drz. fits 


hst_09042_41_wfpc2_f814w_wf_drz.fits 

HST_MOS_1007_ACS_WFC_F814W_drz.fits 
hst _12880_77_wfc3_uvis_f814w_drz.fits 


hst _08602_29 _wfp c2 _f814 w _wf _drz .fits 
hst _08602_30 _wfp c2 _f814 w _wf _drz .fits 
HST_11575_01_ACS_WFC_F814W_drz.fits 


hst _08100_02 _wfp c2 _f814 w _wf _drz .fits 
hst _09042_o 5 _wfp c2 _f814 w _wf _drz .fits 
hst _11570_22_wfc3_uvis_f814w_drz.fits 
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Table A2. Continuation of archival images used 



Date 

Instrument 

Exposure (s) 

At (yr) 

HLAfilename or Dataset 

NGC 4038 

Epoch 1: 

1996-01-20 

WFPC2 

2120 

- 

hst _05962_01 _wfpc2_f814w_wf_drz .fits 

Epoch 2: 

2004-07-21 

ACS/WFC 

1680 

8.51 

HST_10188_10_ACS_WFC_F814W_drz.fits 

Epoch 3: 

2010-01-22 

WFC3 

1032 

14.02 

hst _11577_20 _wfc 3 _u vis _f814 w _drz. fit s 

NGC 4321 

Epoch 1: 

1994-04-23 

WFPC2 

3600 

- 

hst _05397_13_wfpc2_f814w_wf_drz.fits 

Epoch 2: 

1994-06-19 

WFPC2 

3600 

0.16 

hst _05397_12_wfpc2_f814w_wf_drz. fits 

Epoch 3: 

1996-07-29 

WFPC2 

1210 

2.27 

hst _06584_02_wfpc2_f814w_wf_drz.fits 

NGC 4414 

Epoch 1: 

1995-04-06 

WFPC2 

2500 

- 

hst_05397_27_wfpc2_f814w _wf _d r z. fi t s 

Epoch 2: 

1995-06-10 

WFPC2 

2500 

0.18 

hst_05397_2i_wfpc2_f814w _wf _d r z. fi t s 

Epoch 3: 

1999-04-29 

WFPC2 

6400 

4.07 

hlsp_sgal_hst_wfpc2_n4414-l_f814w_vl_mosaic-sci_drz.fits 

NGC 4496 

Epoch 1: 

1994-06-06 

WFPC2 

4000 

- 

hst _05427_02_wfpc2_f814w_wf_drz.fits 

Epoch 2: 

1994-08-06 

WFPC2 

4000 

0.21 

hst _05427_5c_wfpc2_f814w_wf_drz.fits 

Epoch 3: 

2001-03-17 

WFPC2 

640 

6.78 

hst _08599_64_wfpc2_f814w_wf_drz.fits 

NGC 4639 

Epoch 1: 

1995-05-17 

WFPC2 

7800 

- 

hst _05981 _04_wfpc2_f814w_wf_drz .fits 

Epoch 2: 

1995-06-03 

WFPC2 

5200 

0.047 

hst _05981 _08_wfpc2_f814w_wf_drz .fits 

Epoch 3: 

2009-08-07 

WFC3 

800 

14.24 

hst_11570_23_wfc3_uvis_f814w_drz.fits 

NGC 4651 

Epoch 1: 

1994-05-20 

WFPC2 

660 

- 

hst _05375_09 _wfpc2 _f814w_wf _drz. fits 

Epoch 2: 

2007-04-27 

WFPC2 

1200 

12.9 

hst _10803_03_wfpc2_f814w_wf_drz.fits 

Epoch 3: 

2010-11-21 

ACS 

1090 

16.5 

jbks01020_drz.fits 
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